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are determined. 



C^' I. INTRODUCTION 



Singular (or contact or zero-range) interaction potentials were introduced in quantum 
mechanics more than sixty years ago [|1|-|^. Since that time the attitude of physicists and 
mathematicians to this subject was varying, starting from "it is impossible", then to "it 
is evident", and finally arriving at "it is interesting" (for a review see monograph |^). A 
mathematically consistent and rigorous treatment of the subject was developed |^, basing 
on the notion of self-adjoint extension of a Hermitian (symmetric) operator. 

Singular interaction is involved in quantum field theory when, for example, a spinor field 
is quantized in the background of a pointlike magnetic monopole in threedimensional space or 
a pointlike magnetic vortex in twodimensional space. In these cases the Dirac Hamiltonian, 
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in contrast to the Schrodinger one, is free from an explicit 5-function singularity; nonetheless 
the problem of self-adjoint extension of both Dirac and Schrodinger operators arises, albeit 
for different reasons (see, for example ||^). A distinguishing feature is that a solution to the 
Dirac equation, unlike that to the Schrodinger one, cannot obey a condition of regularity 
at the singularity point. It is necessary then to define a boundary condition at this point, 
and the least restrictive, but still physically acceptable, condition is such that guarantees 
self-adjointness of the Dirac Hamiltonian. Thus, effects of polarization of the fermionic 
vacuum in a singular background (such as a pointlike monopole or a pointlike vortex) appear 
to depend on the choice of the boundary condition at the singularity point, and a set of 
permissible boundary conditions is labelled, most generally, by the values of self-adjoint 
extension parameters. In contrast to the Schrodinger case, the extension in the Dirac case 
does not reflect additional types of interaction but represents complementary information 
that must be specified when describing the physical attributes of the already posited singular 
background configuration. 

As a consequence, the fermionic vacuum under the influence of a singular background can 
acquire rather unusual properties: leakage of quantum numbers from the singularity point 
occurs. While in the case of a monopole there is leakage of charge to the vacuum, which 
results in the monopole becoming the dyon violating the Dirac quantization condition and 
CP symmetry |l7|-[TI|, in the case of a vortex (the Ehrenberg-Siday-Aharonov-Bohm potential 
l^ , p!3[| ) the situation is much more complicated, since there is leakage of both charge and 



other quantum numbers to the vacuum. Apparently, this is due to a nontrivial topology of 
the base space in the latter case: tti = in the case of a space with a deleted point, and 
TTi = Z in the case of a space with a deleted line (or a plane with a deleted point); here tti 
is the first homotopy group and Z is a set of integer numbers. For a particular choice of 
the boundary condition at the location of a singular vortex it has been shown that charge 
|jl^,|T5|, current |T6[ and angular momentum |T^ are induced in the vacuum. The induced 



vacuum quantum numbers under general boundary conditions which are compatible with 
self-adjointness have been considered in Refs. ||T8|-pl[]. 



A pointlike static magnetic vortex (the Ehrenberg-Siday-Aharonov-Bohm configuration) 
in 2 + 1 dimensional space-time has the form 
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W(x) = -<l>(°) , 1^2/ N ^(0) 1 , (1.1) 

^^ (Xl)2 + (x2)2' ^""^ "" (2;l)2+(a;2)2' ^^"^^ 

a X V(x) = 27r$(°)(5(x), (1.2) 

where $^°-' is the vortex flux in 27r units, i.e. in the London (27r^ce~^) units, since we use 
conventional units h = c= 1 and coupling constant e is included into vector potential V(x). 
The wave function on the plane (x^, x^) with punctured singular point x^ = x^ = obeys 
the most general condition (see |18| for more details) 

<r,(/? + 27r| =6*2'^^ <r,¥p|, (1.3) 

where r = \J{x^Y + {x'^Y ^^^ ^ — ai'ctan(x^/x^) are the polar coordinates, and T is a 
continuous real parameter which is varied in the range < T < 1. It can be shown (see, for 
example, |1^,|18|) that T as well as ^^^^ is changed under singular gauge transformations. 



whereas difference ^C^) — T remains invariant. Thus, physically sensible quantities are to 
depend on the gauge invariant combination $(°) — T which will be for brevity denoted as 
the reduced vortex flux in the following. 

Thus far the effects of polarization of the massive fermionic vacuum have been studied. 
In the present paper we find the energy density which is induced by a singular vortex in the 
massless fermionic vacuum. In the next section, the most general set of physically acceptable 
boundary conditions at the singularity point x = is defined. In Section III we show how 
the problem of both ultraviolet and infrared divergences in vacuum characteristics is solved 
with the help of zeta function regularization. This allows us to get immediately in Section 
IV the vacuum energy density; also here the effective action and the effective potential are 
considered. We summarize results and discuss their consequences in Section V. The method 
of self-adjoint extension is employed to get the boundary condition at the singularity point 
in Appendix. 



II. QUANTIZATION OF A SPINOR FIELD AND THE BOUNDARY 
CONDITION AT THE LOCATION OF A VORTEX 

The operator of the second-quantized spinor field is presented in the form 

^(x,t) = ^ e"'^'' < x|A > OA + ^ e-*^^* < x|A > 6+, (2.1) 

Ex>0 Ex<0 

where a^ and ax (6^ and bx) are the spinor particle (antiparticle) creation and annihilation 
operators satisfying anticommutation relations 

[ax,a+]+ = [bx,bl,]+=<X\X'>, (2.2) 

and < x|A > is the solution to the stationary Dirac equation 

if <x|A >=Ea <x|A >, (2.3) 

H is the Dirac Hamiltonian, A is the set of parameters (quantum numbers) specifyinff a 

....,...,.„,,.,..„ .,....^......._„..„...... 

integration (with a certain measure) over continuous values of A. The ground state |vac > 
is defined conventionally by equality 

a^jvac >= bx\vac >= 0. (2-4) 

In the case of quantization of a massless spinor field in the background of static vector field 
V(x), the Dirac Hamiltonian takes the form 

H = -ia[d-iV{x.)], (2.5) 

where 

ex = 7°7, P = 7°, (2.6) 

7'^ and 7 are the Dirac 7 matrices. In the 2+1-dimensional space-time (x, t) = {x^,x'^,t) 
the Clifford algebra has two inequivalent irreducible representations which can be differed 
in the following way: 



Choosing the 7*^ matrix in the diagonal form 



7 = 0-3, 



(2. 



one gets 
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,2^ZXsj 



laiC 



-fo-sxs^ 7^ = e 



2'^iXsj 
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0-3 Xs 



(2.9) 



where o"i,o"2 and (T3 are the Pauh matrices, and Xi and X-i are the parameters that are 
varied in the interval < x^ < 27r to go over to the equivalent representations. 

A solution to the Dirac equation (2.3) with Hamiltonian (2.5) in background (1.1), that 
obeys condition (1.3), can be presented as 

^/„(r,E)e^("+^)^ ^ 

\(7„(r,E)e*("+^+^)'^ 



< x.\E,n >-- 



n e z. 



(2.10) 
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where the column of radial functions 



J n 



satisfies the equation 
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e*^= \dr + s{n - $(°) + T + s)r-^] 
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-^^Xs 



-dr + s{n - $(0) + T)r-i] 
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(2.12) 



is the partial Dirac Hamiltonian. When reduced vortex flux $*^°^ — T is integer, the re- 
quirement of square integrability for wave function (2.10) at r ^ provides its regularity, 
rendering the partial Dirac Hamiltonian hn for every value of n to be essentially self-adjoint. 
When $*■''-' — T is fractional, the same is valid only for n 7^ uq, where 



no = [[$(°) -n + l- l-s, 



(2.13) 



[[u]] is the integer part of a quantity u (i.e., the greatest integer that is less than or equal 
to u). For n = Uq, each of the two linearly independent solutions to Eq.(2.11) meets 



the requirement of square integrability at r -^ 0. Any particular solution in this case is 
characterized by at least one (at most both) of the radial functions being divergent as r~^ 
(p < 1) at r ^ 0. If one of the two linearly independent solutions is chosen to have a regular 
upper and an irregular lower component, then the other one has a regular lower and an 
irregular upper component. Therefore, in contrast to operator /i„ (n ^ no), operator /i^^ is 
not essentially self-adjoint [} The Weyl - von Neumann theory of self-adjoint operators (see, 
e.g., Refs. HJ23|]) has to be employed in order to consider the possibility of a self-adjoint 
extension in the case of n = no. It is shown in Appendix that the self-adjoint extension 
exists indeed and is parametrized by one continuous real variable denoted in the following 
by 0. Thus operator hn^ is defined on the domain of functions obeying the condition 

6 TT TT 

cos(s- + -) XmiUrffn^ = -e*^= sin(s- + -) lim(/ir)i-Vo, (2-14) 

2 4 ''^O / 4 r—>0 

where /i > is the parameter of the dimension of inverse length and 

F = .{|$W-T|} + i-i., (2.15) 

{| M 1} = M — [[u]] is the fractional part of a quantity u, < {| m |} < 1; note here that Eq.(2.14) 
implies that < F < 1, since in the case of F = | — |s both /„q and gno obey the condition 
of regularity at r ^ 0. Note also that Eq.(2.14) is periodic in O with period 27r; therefore, 
without a loss of generality, all permissible values of G will be restricted in the following to 
range — tt < O < tt. 

III. ZETA FUNCTION 

In the second-quantized theory the operator of energy is defined as 



^A corollary of the theorem proven in Ref. |22] states that, for the partial Dirac Hamiltonian to 



be essentially self-adjoint, it is necessary and sufficient that a non-square-integrable (at r — > 0) 
solution exist. 



8 = Jd'x^[^+{^,t),H^{^,t)]_ = ^{E.a+a, - E,bp, - ^\E,\), (3.1) 

thus the vacuum expectation value of the energy takes the form 

S =< vac|.^|vac >= — Y^I-EaI = — / d'^xtr < xl \H\ Ix > . (3.2) 

The latter expression is ill-defined due to divergences of various kinds. First, there is a bulk 
divergence resulting from the integration over the infinite twodimensional space. But, even 
if one considers the vacuum energy density, 

^x = -^tr < x| |iJ| |x >, (3.3) 

still it remains to be divergent. There is a divergence at large values of momentum of 
integration, /c — > oo. To tame this divergence, let us introduce the zeta function density 

Cx(-2) =tr<x||i/|-2-|x>, (3.4) 

which is ultraviolet convergent at sufficiently large values of Kez. However, exactly at these 
values of Kez the integral corresponding to Eq.(3.4) is divergent in the infrared region, as 
k ^ 0. To regularize this last divergence, let us introduce fermion mass m, modifying 
definition (3.4): 

C^{z\m) =tr<x\\H\-^'\x>, (3.5) 

where 

H = -ia[d-iy{x)]+pm, (3.6) 

and it is implied that the complete set of solutions to the equation 

^<x|A >=Ea <x|A >, (3.7) 

instead of those to Eq.(2.3), is used. 

In the background of a singular magnetic vortex (1.1) - (1-2) the radial functions of the 
solutions to Eq.(3.7) take the form: 
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In 



v-j 



1 + mE ^Ji^F{kr)e 



«Xs 



l = s{n-no)>0, (3.8) 



y sgn(E) V 1 - mE-^Ji+i^F{,kr) j 



Jn 



( 



2v^ 



1 + mE^^Ji'+Fikr)e 



IXs 



\ 



I' = s{no - n) > 0, (3.9) 



y -sgn(E) V 1 - mE-^Jii_i+F{kr) j 
( -AC, \ 



h 



no 



~(C) 
\ 9na j 



2^7r[l + sin(2z/e) cos(F 



TT 



( -J 



X 



1 + mE 1 [sin(i>^) J_ir(fcr) + cos(z/£.) Ji?(A;r)]e 



«Xs 



\ 



\ sgn(_E)v 1 — m£'^i[sin(z/^) Ji_i7'(A;r) — cos(z/^) J_i+i7'(A;r)] / 



(3.10) 



where A; = v -E^ — m^, Jp(M) is the Bessel function of order p and 



tan(z/^) = sgn(i?)^ 



' 1-mE-i ^ fcy^"' r(l-F) /^ ^ 



(3.11) 



r(M) is the Euler gamma function; note that the radial functions of irregular solution (3.10) 
satisfy condition (2.14) (see Appendix). Note also that Eqs.(3.8) - (3.10) correspond to the 
continuum, \E\ > |m| |. In addition to them, in the case of 



sgn(m) cos 6 < 0, 



(3.12) 



an irregular solution corresponding to the bound state appears. Its radial functions are 



/ ?(BS) \ 
Jno 



~(BS) 
\ 9no I 



sin(F7r) 



/ ^1 + m-^EBsKFiKr)e'^^ \ 



TT V 1 + (2F - l)m-^E_ 



BS 



(3.13) 



y sgn(m) Vl - m-^EBsKi^FiK,r)] J 
where k = ^m^ — E\^^ Kp{w) is the Macdonald function of order p and the bound state 
energy E = Ebs (I-^^bsI < |"^|) is determined implicitly by the equation 

(l + m-i^Bs)^"^ . . /|m|\^^~^r(l-F) , /e 7i' 



l-m-i^RsF 



|m|\'^~'r(l-F) 

-sgnlml I ; — - — tan s 

^ ^ '2^7 r(F) V 2 4 



(3.14) 



^In a 2+1-, as well as in any odd-, dimensional space-time mass parameter m in Eq.(3.6) can take 
both positive and negative values. 
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Regular solutions (3.8) and (3.9) yield the following contribution to zeta function density 
(3.5): 



oo 



[Cx(^|m)]REG = ^J dkk\E\-'^ J2 {J2l^l + mE-')Jlp{kr) + {1 - mE-')J^^,_p{kr)] + 

sgn{E) ^=1 

CO 

+ J][(l + mE'')jf,^p{kr) + (1 - mE'')Ji,_i+F{kr)]}. (3.15) 

i'=i 

Summing over the energy sign and over I and /', we get the expression 

oo kr 

[Cx(^|m)]REG = - [dkk\E\-'^ f^[F4{y) + {l-F)JlM], (3.16) 

TT y J y 





which in the case of Re z > 1 is reduced to the form 



1 r (jk ~ 

[Cx(^|m)]REG = ^^(^_^^ j ^\Er'^{FJl{kr) + (1 - F)Jl_,{kr)]. (3.17) 



Irregular solution (4.10) yields the following contribution to Eq.(3.5): 

oo 

[Cx(^|m)]iRREG = ^ J dkk\E\-'''^{Af,'-'^e^m+) - L^^)]J'_^{kr) + 



+A^'~^^k~^^'~^\im-\E\fL^+) - (m+ |E|)2L(_)]Ji%(A;r) + 
+2[(m + |^|)L(+) - (m - |^|)L(_)] J_^(A;r) J^(A;r) + 2[(m - |^|)L(+)- 
-(m + |^|)L(_)] Ji_^(A;r) J_i+^(A;r) + A-'fi^^-'k~^^[{m + |^|)'L(+)- 



m 



|^|)2L(_)] j|(A;r) + A-V'''-'A;2^'-^H^(+) " ^(-)]^!i+F(A:r)}, (3.18) 

where summation over the energy sign has been performed and 

L(±) = [A/-'^A;-2(i-^)(_rri ± |E|) + 2 cos(F7r) + A-'fi^^-'k-^^im ± 1^1)]"^ (3.20) 
The contribution of bound state solution (3.13) to Eq.(3.5) is the following: 

[Cx(;^|m)]Bs = ^^^ fl^''fl7 n [(^ + EBs)KUKr) + {m- E^s)KlA'^r)]. (3.21) 
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By deforming the contour of integration in the complex w-plane, Eq.(3.17) in the case 
of 1 < Re z < 2 is transformed to the following expression 

r> / I M ImP*^^"^) sin(2;7r) r,,^ .-. 

2tT{Z — 1) 7T^{Z — 1) 

00 

X [ —(w^ - m^r^y-'[FlF(w)KF(w) + (1 - F)h_p(w)Ki_F(w)l (3.22) 

J w 

\m\r 

while Eq.(3.18) in the case of | < Rez < 1 is transformed to the following one 

00 
[Cx(^|m)]iRREG = ^^^r2(-i) j dww{w' - m'r'r'[lF{w)KF{w) + hMw)Ki-FH] + 

\m\r 

2sin(F7r) . , , 2(z-^) 
X / dww[w — m r ) - 



\m\r 

^ ^ + ^s(2F - 1) ^^^ + EBs)KU^r) + {m- EBs)^?_^(«:r)]; (3.23) 

here Ip{w) is the modified Bessel function of order p. The integral in Eq.(3.22) can be 
analytically continued to domain | < Rez < 2. In the case of | < Re^; < 1 this integral is 
decomposed into two terms: 

|^|2{l-.) 



[Cx(^|"^)] 



REG 



27r{z - 1) 
00 



\m\r 

00 

I 2sin(F7r) sin(^^2(.-i) 

TT^ Z — 1 



00 

f dw{w^ - m^r'^f-'KF{w)Ki_F{w), (3.24) 



\m\r 

the last of which can be analytically continued to domain Re 2; < 2. Note also that the 
second integral in Eq.(3.23) can be analytically continued to domain Re^; < 1. 
Summing Eqs.(3.21), (3.23) and (3.24), we get 



A / I X |mP(^ ''^ 2sin(F7r) sin(2;7r) ■,,-,,. , , , ., ., .,,-, ,^^,.^^ , , 

C^{z\m) = ^ + ^ — '- — ^Z^A^-i) / dw{w' - m'r'y-'KF{w)Ki^F{w) + 

27r[Z — 1) TT'^ z — I J 

\m\r 
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2 sin(F7r) 



sm 



VT" 



{zny^'-'^x 



X / aww{w — m r ) - 



w) 



(3.25) 



|m|r 



i.e., the terms which are defined only in domain | < Rez < 1 are cancelled. 

Note that the first term in Eq.(3.25) is identified with the zeta function density in the 
noninteracting theory (i.e. in the absence of any boundary condition and any background 
field): 



Cf(z|m) 



771 



2(1-2) 



(3.26) 

In the noninteracting theory all vacuum values are simply omitted due to the prescription 
of normal ordering of the product of operators (see, for example, 1^3 )• Therefore, one has 
to subtract Q from (^x ioi the reasons of consistency. Doing this and removing the infrared 
regulator mass, we obtain the renormalized zeta function density 





)-e(^m)] 


= 


^-'-'^-1^14 


r F(i-F)] 

12^ ' 
\-z 


r(F-^)r( 


+ / dww^''^'' [-^l(w) - -^i-f(^)] tai 


ah 


(2F - 1) In 


— +lnA 



(3.27) 



IV. ENERGY DENSITY AND EFFECTIVE POTENTIAL 



Recalling the formal expressions for the vacuum energy and zeta function densities, 
Eqs.(3.3) and (3.4), one can easily deduce that the physical (renormalized) vacuum energy 
density is expressed through the renormalized zeta function density at z - 



1. 

'2' 



s: 



_ /-ren/ _^ 

-2^ I 2' 



(4.1) 



In the background of a singular magnetic vortex (1.1) - (1-2), using Eq.(3.27), we get the 
expression 



SI 



sin(F7r) 



i-F 



2TTr^ 1 6 cos(F7r) 



F(l-F) 



+ 
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oo 

+^ f dww'^\Kl{w) - Ki%(«;)]tanh[(2F - 1) ln(— ) + In A]|. 
'R J fir J 



(4.2) 



At noninteger values of reduced vortex flux $*^°^ — T (i.e. at < F < 1) vacuum energy 
density (4.2) is positive. At half-integer values of the reduced vortex flux (F = ^) we get 

1 



In the case of cos 9 = we get 



cren 



tanCFV) 1. 

47rr3 ^ 2' 



creni 

^- 1^=1 ~ 247r2r3' 



-F(l-F)--T-iF--) 
3 ^ ^ 4 2^ 2^ 



(4.3) 



e = ±s 



TT 



(4.4) 



If cos O 7^ 0, then at large distances from the vortex we get 

tan(F7r 



£1 



Airr^ 



<F--) 



iF(l-F)-i + V-i| 

3 ^ ^4 2' 2' 



(4.5) 



Going over to imaginary time t = —ir, let us consider the effective action in 2+1- 
dimensional Euclidean space-time: 



5(f+i)[V(x)] = -ln{A^-^ / c^^c^^+exp[- / dTd'^x'^+{-i(3dr ~ i(3H)^} 



\nl)et[{-il3dr - if3H)m 



™-ii 



(4.6) 



here A^ is a normalization factor, parameter m is inserted just for the dimension reasons, 
while fermion mass m (see Eq.(3.6)) is introduced in order to tame the infrared divergence. 
The real part of the effective action is presented in the form 

1 



Re5^2%[V(x) 



drd X tr < x, t\ \n[{-d^ + H )m ] |x, r > 



(4.7) 



Let us define the zeta function density in threedimensional space {x\x ,t): 



Cx,r(-2|m) = tr < X, t\ {-d^ + if^)"^|x, r > . 



(4.8) 



^The imaginary part of the effective action vanishes in the case of a static background. 
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Then an ultraviolet regularization of Eq.(4.7) can be achieved by expressing its integrand 
through Eq.(4.8): 

--tr < X, r I ln[(-92 + H')m-^] |x, r >= - [—(^^,{z\m)]\^^, + -Cx,r(0|m) In ml (4.9) 

In the background of a singular magnetic vortex (1.1) - (1-2) we get, similarly to 
Eq.(3.25), the following expression 



sin(F7r) 

7 
712 



cosUtt 



Cx,t(^|"^) = 

r(^) 

sin(F7r) 



\m 



3-2z 



Tiz 
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V{z\ 



\m\r 



Ajj, 



1-2F ( w \2F 1^2 



Kl{w) + A-^ii^^-\^f^^~^)Kl_p{w) 



(4.10) 



note that the first term in Eq.(4.10) corresponds to the case of the noninteracting theory: 



C(^l^) 



772 



3-2z 



rf^- 



A-K2 T{z) 



(4.11) 



Note also relation 



Cx,r(0|m) = 0, (4.12) 

which ensures the independence of the effective action on m^; thus Eq.(4.9) takes the form 



\ra\r 



1 

-1 

2 



\m\r 

The effective potential in the massless theory is defined as 



. (4.13) 



1 

2 m—*0 



W^U, t) = -^ lim tr(x, t\ In 






|x,r). 



(4.14) 
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Defining the renornialized zeta function density 



Ko), 



C(^) = lim [Cx,.(^|m) - a'>iz\m)], (4.15) 

we get, similarly to Eq.(4.9), 

^^''(x, r) = ^[^C^(^)] |.=o +lc:m Inm^ (4.16) 

In the background of a singular magnetic vortex (1.1) - (1-2), using Eq.(4.10), we get 



^^^„, , sin(F7r) , T{z—-::) „ „ 

2tt2 i [z) 



nVa-z) 



\ 4F(1 - F) 
3-2z 



Ti^-z + F)T{^-z-F) + 





in particular 



4 r(2-;2) 

oo 

+ f dww^^^-'\KUw) - K^_piw)]tanh[{2F - I) \n{—) + In A] I , (4.17) 

J l^'^ I 



C.^(0) = 0, (4.18) 



and, thence, we arrive at the remarkable relation 

W\-K,T)=8r. (4.19) 

where Ell'^ is given by Eq.(4.2). 

Although the last relation looks rather natural and even evident, let us emphasize here 
that it is a consequence of the relation between the renornialized zeta function densities of 
different spatial dimensions. 



[^o^)]L=o = -cr(-^ 



0^)]L=o = -Cr(-7T), (4-20) 



and relation (4.18). As it has been shown in Ref. [^, relation (4.20) can be in general broken 
in spaces of higher dimensions. Moreover, both left- and right-hand sides of Eq.(4.20) can 
have nothing to do with the true vacuum energy density. Fortunately, this is not relevant 
for the case considered in the present paper, and, indeed, in the background of a singular 
magnetic vortex the vacuum energy density coincides with the effective potential. 
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V. CONCLUSION 

In the present paper we show that the massless fermionic vacuum under the influence of 
a singular magnetic vortex (1.1) - (1-2) in 2+1-dimensional space-time attains the energy 
density (4.2) which decreases at large distances from the vortex as inverse power with integer 
exponent (4.5). 

The most general set of boundary conditions at the location of the vortex is used, see 
Eq.(2.14), providing the self-adjointness of the Dirac Hamiltonian; thus the vacuum energy 
density is depending on self-adjoint extension parameter 6 or A (3.19). As to the depen- 
dence on vortex flux $(°\ it has been already anticipated in Introduction that all vacuum 
polarization effects are gauge invariant and thus depend on reduced vortex flux $(°'' — T 
rather than on $('^^ or T separately. Note also that at half-integer values of $*^°-' — T (i.e. at 
F = i) the vacuum energy density is given by Eq.(4.3). 

Among the whole variety of boundary conditions which are specified by self-adjoint ex- 
tension parameter 9, condition cos 6 = (or 6 = ±|) is distinguished, since it corresponds 
to one of the two components of a solution to the Dirac equation being regular for all n: 
if 6 = s|, then the lower components are regular, and, if 6 = — s|, then the upper com- 
ponents are regular. This condition is parity invariant, and under it the vacuum energy 
density is given by Eq.(4.4). It should be noted that this condition is extensively discussed 
in the literature, being involved into the two most popular ones: the condition of maximal 
simplicity [^] 



sf, s(<l>(°)-T)>0 

e = { ) (5.1) 

-sf, s(<l>(o)-T) <0 



and the condition of minimal irregularity [p!5|,pi 
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0= < 



s ^ , 



0, 



-l<s({|$W-T|}-l)<0 
{|$(o)_T|} = i 



> ; 



-sf, 0<s({|'^^°^-T|} 



-) < - 



(5.2) 



here, both in Eqs.(5.1) and (5.2), it is imphed that {|$(°) - T |} 7^ 0. 
Under condition (5.1) we get 



pren 



Wflg!Lz2nh) {| $(0) _ T |}(i _ {| ${o) _ T |}2), $(0) _ T > 

W(l^l!lz]n}M (1 _ {| $(0) _ T |})[i - (1 - {| $(0) - T \})% $(0) - T < 



(5.3) 



Under condition (5.2) we get 

' W{|!i^]qM({| $(0) _ T 1} - i)[3|{| $(0) - T 1} 



s: 



-2({|$(o)-T|}-i)2-l], 



1 

247r2r3 ' 



{|$W-T|}^i 



(5.4) 



It is clear that Eq.(5.4), in contrast to Eq.(5.3), is periodic in the value of the vortex 
flux. 

As it should be expected, the vacuum energy density is invariant under transitions to 
equivalent representations of the Clifford algebra (i.e. independent of Xs)- H should be 
emphasized that the vacuum energy density is also invariant under the transition to an 
inequivalent representation (i.e. under s -^ — s). 
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APPENDIX 

Let us consider a general case of massive Hamiltonian H (3.6). The relevant partial 
Hamiltonian has the form 



"no 



( 



\ 



e-'^H-dr - Fr-^) -m 



J 



(A.l) 



Let h be the operator in the form of Eq.(A.l), which acts on the domain of functions ^^{r) 
that are regular at r = 0. Then its adjoint /i^ which is defined by the relation 



drr[h^^{r)]^^\r) 



drr[^{r)y[hC{r)] 



(A.2) 



'0 Jo 

acts on the domain of functions ^(r) that are not necessarily regular at r = 0. So the 
question is, whether the domain of definition of h can be extended, resulting in both the 
operator and its adjoint being defined on the same domain of functions? To answer this, 
one has to construct the eigenspaces of h'^ with complex eigenvalues. They are spanned 
by the linearly independent square-integrable solutions correspoding to the pair of purely 
imaginary eigenvalues. 



h^^^{r) = ±iixe{r), 



(A.3) 



where /i > is inserted for the dimension reasons. It can be shown that in the case of 
Eq.(A.l) only one pair of such solutions exists, thus the deficiency index of h is equal to 
(1,1). This pair is given by the following expression 



^^(r) = — 



V 



e*^= exp [±^sgn(m)ri]K p(fir) 
sgn(m) exp [^^sgn{m)ri]Ki_F{fir) 



(A.4) 



where A^ is a certain normalization factor and 



u = V /i^ + m'^, ri = arctanf - — r ] . 

\\m\/ 



(A.5) 



Self-adjoint extended operator h^" is defined on the domain of functions of the form 
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J no 
\y9no j 



e° + c(e+ + e^V), 



(A.6) 



where c is a complex parameter and Qg is a real continuous parameter which depends, in 
general, on the choice between the two inequivalent representations of the Clifford algebra. 
Using the asymptotics of the Macdonald function at small values of the variable, we get 



or 



is \ 

J no 






e'^' cos {\[0s- sgn(m)r7]} 2^r(F)(/ir)-^ 

sgn(m) cos {\[0s + sgn(m)?7]} 2^^^r(l — F)(/ir)~-'^+-^ 



tan \-6s sg\i{m)ri] sin?] — sgn(m) cos?] > lim(/ir) /, 



r^O 



/no 



(A.7) 



-^"^2--^;^lim(;ir)--^„„. 



(Ai 



Defining new parameter O by means of relation 



, TT 

tan s 1 — 

2 4 



11 1 " /2u\2^-i ViF) 

tan[-6'^ - -sgn(??2)?]]sin?]-sgn(??2)cos?]^ \^-rj . . , (A.9) 



we get Eq.(2.14). 

Certainly, both 9s and O can be regarded as self-adjoint extension parameters which are 
to specify the boundary condition at r = 0. The use of O in this aspect may seem to be 
more preferable just for the convenience reasons, because Eq.(2.14) looks much simpler than 
Eq.(A.8). In particular, Eq.(2.14), in contrast to Eq.(A.8), is independent of ??2 and remains 
explicitly invariant under s -^ —s (0 is independent of s). 

In the limit of ??2 -^ ±0 Eq.(A.9) takes the form 

TT^ 



tan s 



4 



tan [Us + ^sgn(?r?)7r]22^-i^;^^^ 



(A.IO) 



then we get 



e., = s9. 



(A.ll) 



where 6 is independent of s. 

Concluding this appendix, let us note that all characteristics of the massless fermionic 
vacuum are depending on A (3.19) rather than on itself, and A is expressed through 9 in 
the following way 

A = -tan [-s6l + -sgn(m)7r]. (A.12) 
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